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J> ■ The Antarctic Muon And Neutrino Detector Array (AMANDA) has been tak- 

Q\ ' ing data since 2000 and its data acquisition system was upgraded in January 

2003 to read out the complete digitized waveforms from the buried Photo- 
, multipliers (PMTs) using Transient Waveform Recorders (TWR). This system 

currently runs in parallel with the standard AMANDA data acquisition sys- 
f^*^ ' tern. Once AMANDA is incorporated into the 1 km 3 detector IceCube, only the 

i TWR system will be kept. We report results from a first atmospheric neutrino 

f**^ , analysis on data collected in 2003 with TWR. Good agreement in event rate 

and angular distribution verify the performance of the TWR system. A search 
' of the northern hemisphere for localized event clusters shows no statistically 

significant excess, thus a flux limit is calculated, which is in full agreement with 
^ , previous results based on the standard AMANDA data acquisition system. We 

» I ■ also update the status of a search for diffusely distributed neutrinos with ultra 



high energy (UHE) using data collected by the TWR system. 
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H , Introduction 



Neutrinos are the only high energy particles able to propagate undeflected 
and unattenuated from the furthest reaches of the Universe. Extragalac- 
tic UHE 7-ray astronomy falters for energies greater than a few tens of 
TeV due to interactions with infrared and Cosmic Microwave Background 
photons. The information carried by the neutrino messengers from dis- 
tant, unexplored regions of the universe may help to unravel longstanding 
mysteries associated with the origin of the highest energy cosmic rays. Sev- 
eral models^' ^' ^' • ^' ® predict high energy charged particle and neutrino 
emission; the diffuse z^-flux predictions by these models are constrained by 



*http:/ /www. icecube. wisc.edu 
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the observed cosmic ray fluxes'''. AMANDA^, the first neutrino telescope 
constructed in transparent ice, is deployed between 1500 m and 2000 m 
beneath the surface at the geographic South Pole. It is designed to search 
for neutrinos that originate in the most violent phenomena in the observ- 
able universe. Galactic objects like Supernova Remnants (SNR) and extra- 
galactic objects such as Active Galactic Nuclei (AGN) are expected to be 
the powerful engines accelerating protons and nuclei to the highest ener- 
gies, which eventually interact to generate neutrinos. We have searched for 
point sources in the northern sky and for a diffuse flux of UHE (E„ > 10 15 
eV) neutrinos of cosmic origin. We have performed the analysis using data 
collected in 2003 with a new data acquisition system based on Transient 
Waveform Recorders (TWR-DAQ), which we compared to data collected by 
the standard AMANDA readout system (/x-DAQ). 




Fig. 1. The upgraded AMANDA data acquisition electronics with Transient Waveform 
Recorders, here displayed in a system of 2 racks with 72 TWR modules. 



1. The TWR System 

The data acquisition electronics of the AMANDA-II detector (Fig. 1) was 
upgraded in 2003 to read out the complete digitized waveform of the pho- 
tomultiplier tubes (PMTs) using Transient Waveform Recorders (TWR)^. 
The transition was made to run the fi-DkQ}® and TWR-DAQ in parallel to 
verify that the performance of the latter is as good or better. To compare 
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the two systems, the data from 2003 has been analyzed with data sets from 
both readout systems. The performance of the TWR-DAQ is verified by com- 
paring the results for the absolute rate of atmospheric neutrinos and the 
cos{9) distribution with the results from the standard fi-DKQ^. Extending 
the analysis tools to include TWR data required several new developments: 
(1) The TWR-DAQ measures the integrated charge Q from the full waveform 
of the PMT pulses, in contrast, the /z-DAQ only measures the maximum 
amplitude of the PMT pulse in a 2 /is window. (2) Various time offsets 
were taken into account for the TWR-system, none of which impacted the 
/i-DAQ. A time resolution of few nanoseconds was extracted. (3) The size 
of the TWR-DAQ data set (20TB for 2003) by far exceeds that of the fi-DkQ 
(~1TB), requiring an improved data handling tools. 

2. Analysis 

The data information of the two systems has been merged according to 
GPS time and the fraction of overlapping PMTs participating in the event 
in both systems. For the point source analysis we restricted the capabilities 
of the TWR-DAQ system to mimic the features of the /x-DAQ system as closely 
as possible. The timing and amplitude information extracted from wave- 
forms has been used as input parameters to perform PMT-pulse cleaning, 
TOT (Time-Over-Threshold) and crosstalk a cleaning. A Gaussian fit of the 




Fig. 2. (Left) Time difference between the PMT pulses recorded by the two acquisition 
systems At = t M — txwR- (Right) Calibrated amplitude normalized to 1 photo-electron 
(p.e.) value for /i-DAQ and TWR-DAQ data. See text for details. 



a Any phenomenon by which a signal transmitted on one channel of a transmission system 
creates an undesired effect in another channel. 
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distribution for At = — txwR yields a At — 4.2 ns (Fig. 2 (left)), which is 
dominated by the systematic error of the time jitter between independent 
flash ADC clocks of the TWR system. The TWR-DAQ timing calculations are 
relative to the values measured by the /i-DAQ. The timing of the TWR-DAQ 
system includes two sources of jitter, the digitization window of 10 ns and 
the time fluctuation of the flash ADC clocks within the same 10 ns interval. 
This accounts for the a At ~ 4 ns b . Amplitudes are also calibrated by ex- 
tracting the number of photo-electrons (Npe) detected from the peak ADC 
of the /i-DAQ and charge Q of the TWR system and normalized to 1-pe am- 
plitude. By integrating the charge from pulses in the waveform, the dynamic 
range of the TWR-DAQ extends to Npe of ~ 100, a factor of three higher than 
the /z-DAQ. Fig. 2 (right) shows the reconstructed amplitude of the TWR-DAQ 
compared to the /i-DAQ, which indicates a stable power law distribution 
extending up to 100 Npe, while the /i-DAQ system shows a "knee" around 
10 Npe. The knee is due to the amplitude saturation of channels read out 
by optical fibers, which comprise ^40% of the AMANDA channels. Af- 
ter cleaning, the muon track is reconstructed from the pulse times of the 

1 

PMTs. Details on the reconstruction techniques can be found in . Table 1 



Table 1. Passing rates for increasing level of 
data selection criteria for the TWR-DAQ and 
fi-DAQ data analysis. 



Selection 


TWR-DAQ 


/i-DAQ 


Lcvcl-0 (Raw) 


1.86 x 10 9 


1.86 x 10 9 


Level- 1 


1.25 x 10 8 


1.25 X 10 8 


Level-2 


2.56 x 10 6 


1.99 X 10 fi 


Level-3 (Final) 


1112 


1026 



summarizes the passing rates from the raw data level to the final sample 
of atmospheric neutrinos. The event selection criteria used in this analysis 
follows the method described in 13 . Fig. 3 shows that the angular mismatch 
for the final neutrino sample gaq is 0.5° where A6 = 9^ — Otwr- This value 
is expected from studies of the precision of the global minimizcr in the re- 
construction program, by reconstructing the same event sample twice with 
a 32-iteration algorithm. The azimuth angular mismatch <ja4> is 0.7°, which 



b For year 2005 the phases of the TWR modules were synchronized to avoid the time 
jitter between independent flash ADC clocks. 

c High voltage values were lowered in January 2005 to increase linear dynamic range of 
optical channels. 
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Fig. 3. (Left) The zenith AS = M — 9twr difference distribution between the fi-DkQ 
and the TWR-DAQ systems, (right) the azimuthal A<j> = 4>ti — 4>TWR difference distribution. 



3. Search for Point Sources of Neutrinos 

From the analysis based on TWR-DAQ data, 1112 neutrinos are observed com- 
pared to 1026 neutrinos from the /i-DAQ data analysis. The small differences 
in the event rate are compatible with the small differences in analysis pro- 
cedures described in Section 2. Fig. 4 shows the cos(0) distribution of the 
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Fig. 4. Distribution of cos(6>) after final selection representing the atmospheric neutrino 
sample observed from the TWR-DAQ and the /i-DAQ analyses. 
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Fig. 5. Sky plot displayed in equatorial coordinates, Right Ascension (RA) and Dec- 
lination (DEC), of the final sample of 1112 atmospheric neutrinos observed from the 
TWR-DAQ analysis. 



atmospheric neutrino sample extracted from the TWR-DAQ and ^t-DAQ data 
analysis. Satisfactory agreement can be seen for the cos(#) distribution 
of the atmospheric neutrinos samples obtained by the different analyses. 
Fig. 5 displays the sky map plotted in equatorial coordinates of the 1112 
atmospheric neutrino candidates observed from the TWR-DAQ analysis. All 
events are distributed in the norther hemisphere with very few events re- 
constructed at the horizon. In order to distinguish if the observed events 
follow a random distribution as expected by atmospheric neutrinos, or are 
an indication of localized event cluster as expected by a source of neutrinos 




Fig. 6. Sky map plotted in coordinates of right ascension and declination of the 1112 
atmospheric neutrino candidates from the TWR-DAQ analysis. The scale on the right reflects 
excess or deficit in terms of the positive or negative deviation with respect to the mean 
background events. 
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of extraterrestrial origin, an analysis is required to estimate the statistical 
significance of the observed events. A full search of the northern sky was 
performed to look for any localized event cluster. The full scan is extended 
to 85° in declination, since the limited statistics in the polar bin prevents 
an accurate estimate of the background. Fig. 6 shows the sky map of the 
calculated significance for all observed 1112 events in terms of the a pos- 
itive or negative deviation with respect to the mean background events. 
All observed regions with the highest significance are compatible with the 
background hypothesis. The highest significance observed shows a positive 
deviation of 4.3 a, which corresponds to a probability of 23% for a search 
performed on events with randomized right ascension. 




Fig. 7. (Left) Muon neutrino effective area as a function of E„ M for point source searches 
displayed for different declinations 6. (Right) Average upper limit as a function of sin(f5) 
on Ufj, for a E -2 neutrino flux spectrum from the TWR-DAQ and the /i-DAQ analyses. 



Numerous studies have been performed to implement the best possible de- 
scription of the TWR system in the Monte Carlo simulation. In particular, 
a correct description of the waveform hardware response is now available 
and proper TDC time windows as well as amplitude calculations have been 
now implemented. Therefore a calculation of the neutrino effective area for 
point source search using TWR data is now feasible. Fig. 7 (left) shows 
the effective muon neutrino area as a function of E u computed at the fi- 
nal selection of this analysis for four separated ranges of the declination. 
The Ag^-j extends over six order of magnitudes as ^-energy increases, while 
A-eff decreases as declination increases due to the neutrino absorption in 
the Earth. The average upper limit on as a function of declination is 
shown in the right panel of Fig. 7. The average neutrino flux upper limit 
is determined from the ratio of the average Feldman and Cousins^ upper 
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limit < fi 90 > computed according to the expected mean background and 
observed events, and the number of the expected signal events for a neu- 
trino flux E^d^/dE^ = GeV s _1 cm~ 2 from a point source at declination 
S. Together with the average upper limit extracted from the TWR-DAQ data, 
the limit from the /i-DAQ data of year 2003 is also represented for com- 
parison. These results are comparable, however the TWR-DAQ limit is ~ 10% 
worse than the ^t-DAQ limit, because conservative event selection criteria 
have been applied that do not use the full waveform information in an op- 
timal way. The preliminary average upper limit from the TWR-DAQ analysis 
on the muon neutrino flux with spectrum d<I>„/dE„ oc E~ 2 is E^d^/dEj, 
< 1.8 x 10~ 7 GeV s _1 cm -2 at 90% confidence level, in the energy range 
1.26 TeV < E„ < 1.6 PeV. This upper limit does not include systematic 
uncertainties. AMANDA has submitted a publication on a 5-year search 
for point sources (2000-2004) based on the /i-DAQ data^, which improves 
the limit by ^67% compared to the results from this analysis. 




Fig. 8. (Left) one of the 1112 atmospheric neutrino event detected in the experimental 
array. (Right) a high-energy MC simulated event. The size of the circles represent energy 
deposited in the detector, the colors represent the time profile of photons propagating 
in the ice which have scattered to the PMTs. 
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4. Search for UHE Neutrinos 

The search for UHE neutrinos differs from the point source mainly in two 
aspects. First, the search for point sources is restricted to upward going 
events of energy within the TeV scale, the UHE analysis extends the search 
beyond PeV energies. Second, while point sources are localized in a small re- 
gion of the sky, UHE diffuse neutrinos are expected to come from the entire 
visible sky (2ir sr), with a higher probability to originate from the horizon. 
Fig. 8 shows one of the 1112 upward going neutrino events and an almost 
horizontal UHE simulated neutrino event with energy E ~ 3 x 10 15 eV. 
Simulated events with energy above one PeV deposit a substantial amount 
of light which is recorded by almost all PMTs. In order to better separate 
background events from signal events, we developed new variables which 
exploit the information of the full waveforms. The capability of improved 




Fig. 9. Contour plots displaying the capability of correctly reconstructing the total 
number of photon-electron injected into the detector (left) for the TWR-DAQ and (right) 
for the fi-DAQ system. 



photo-electron reconstruction of the TWR-DAQ system can be seen in Fig. 9, 
which displays a contour plot in logio-scale of the sum of the reconstructed 
number of PE, NPE reco versus the NPE true generated by the signal MC 
simulation. While the /x-DAQ system quickly saturates on the order of few 
thousands of PE's, the TWR-DAQ, by using the after-pulse information, can 
see the millions of PE's which are typical for UHE events triggered by sin- 
gle muons. Fig. 10 summarizes the results of several AMANDA searches 
for diffuse neutrinos at different energy ranges using yU-DAQ data. The ex- 
perimental limits assume a 1:1:1 ratio of neutrino flavors at the Earth due 
to oscillation. The dotted and dashed lines represent a sample of model 
predictions^' ^ : ^' ^ : adjusted for oscillation if necessary. The horizontal 
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solid lines represent limits on the integrated neutrino flux from the diffuse 
analyses with AMANDA-Bf 15 , with AMANDA-II 16 , from the cascade 

17 id in 

analysis , and from the UHE analysis ' , with the NT200 neutrino 
telescope 21 ^ , and at the highest energies from ANITA-lite 21 and RICE 22 . 
These are the most stringent flux limits at 90% C.L. for an E -2 spectrum 
to date. The AMANDA limits have been determined by analysing data 
collected with the /i-DAQ system. However, as it has been shown above, by 
performing a new analysis approach with TWR-DAQ data, it is possible to 
develop new techniques and to improve the current experimental limits for 
the UHE neutrino search. 




4 5 6 7 8 9 10 11 12 13 14 



log (E v [GeV]) 

Fig. 10. Experimental limits on the integrated neutrino flux and detector sensitivity 
for a neutrino flux of E -2 spectrum. The curves represents predictions from theoretical 
models, the solid lines show current AMANDA (all-flavor limits), NT200, ANITA and 
RICE experimental limits. 



5. Conclusion 

The point source search provides the first detailed evaluation of the perfor- 
mance of the TWR-DAQ system. This analysis demonstrated that the TWR-DAQ 
produces similar event rates and angular distribution as the data collected 
by the standard ^-DAQ system. No statistically significant excess has been 
observed after performing a full search of the northern hemisphere of the 
sky for localized event clusters, therefore a flux limit based on TWR-DAQ 
data is calculated. The AMANDA TWR readout is now being incorporated 
into the IceCube data acquisition system. By exploiting the information of 
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the full waveforms from the TWR-DAQ system, it is possible to develop new 
analysis techniques for an improved search for diffuse UHE neutrinos. Cur- 
rently the AMANDA experiment has placed the most stringent neutrino 
flux limits to date, which can be further improved by analyses performed 
with TWR-DAQ data. 
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